Evolution of /nl to the adiabatic limit 
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Abstract. We study inflationary perturbations in multiple-field models, for which £ typically 
evolves until all isocurvature modes decay — the "adiabatic limit" . We use numerical methods 
to explore the sensitivity of the nonlinear parameter /nl to the process by which this limit is 
achieved, finding an appreciable dependence on model-specific data such as the time at which 
slow-roll breaks down or the timescale of reheating. In models with a sum-separable potential 
where the isocurvature modes decay before the end of the slow-roll phase we give an analytic 
criterion for the asymptotic value of /nl to be large. Other examples can be constructed 
using a waterfall field to terminate inflation while /nl is transiently large, caused by descent 
from a ridge or convergence into a valley. We show that these two types of evolution are 
distinguished by the sign of the bispectrum, and give approximate expressions for the peak 
/nl- 



1 Introduction 



In single-field inflation with canonical kinetic terms, the curvature perturbation produced 
at horizon crossing is conserved with nearly Gaussian statistics [1, 2]. Multiple- field models 
support a richer phenomenology, driven by a flow of power from isocurvature modes into the 
curvature perturbation. This flow is sourced dynamically, and where only canonical kinetic 
terms are present the dynamics are determined by the potential. Therefore, there is some 
hope that we may use one to learn about the other. 

How much information could be extracted? A general potential is a complicated land- 
scape, and it is unlikely that observations will be sufficient to single out a specific shape. But 
by piecing together clues from dynamical evolution it may be possible to obtain information 
about the topography of the landscape in our local neighbourhood. This is a form of potential 
reconstruction [3-5]. 

Sensitivity to dynamical effects is helpful when distinguishing observational outcomes. 
Unfortunately, it complicates the task of extracting predictions. In principle, the statistics 
of the curvature perturbation should be tracked until the time of last scattering — where the 
microwave background anisotropy was imprinted — and in our present state of ignorance this 
is an impossible undertaking. Therefore, to connect the physics of inflation with observation, 
we must rely on conservation: if the isocurvature modes are exhausted, quenching the flow 
of power into the curvature perturbation, it will cease to evolve. It is the statistics which 
apply at the onset of conservation which will be inherited by observable quantities. 

This point of view was developed soon after multiple-field models entered the litera- 
ture [6]. For practical purposes we require a characterization of the conditions under which £ 
becomes constant. In the absence of isocurvature modes, conservation of £ was demonstrated 
by Rigopoulos & Shellard [7], Lyth, Malik & Sasaki [1] and later Langlois & Vernizzi [8-12] 
using a gradient expansion. More recently, Naruko & Sasaki applied these arguments to 
noncanonical models [13], where conservation can be subtle [14, 15]. Weinberg developed 
a different approach [16-19], adapting the techniques of Goldstone's theorem to show that 
£ would become massless on superhorizon scales, admitting a time-independent solution. 
Whether this solution is selected is a model-dependent question. 

The statistics of £ are fossilized in the radiation fluid at last scattering, and its two-point 
correlations have been studied since their presence was confirmed by COBE [20, 21]. More 
recently, sophisticated Cosmic Microwave Background (cmb) experiments have raised the 
possibility that three- and higher n-point correlations may be detectable [22, 23]. These cor- 
relations are interesting because, in principle, they are sensitive probes of physical processes 
and dynamical conditions in the early universe [2, 24-27]. But precisely because of these 
desirable properties, such observables carry an unavoidable risk: they may be equally sensi- 
tive to subsequent dynamics, including the process by which £ becomes conserved. We may 
learn important physics from studying the details of this process — but it need have little to 
do with our theories of the very early universe, and if we wish to use n-point correlations to 
study these theories then we should proceed with caution. For this reason it is important to 
understand which predictions of early-universe models can be connected reliably to late-time 
observations. 
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Our imprecise knowledge of physics above the TeV frontier means it is not possible to 
give a complete answer. In this paper we pursue a more modest objective. Focusing on three- 
point correlations — with amplitude measured by /nl [22] — and considering models where the 
flow of power from isocurvature modes is quenched at or near the end of inflation, we study 
how the asymptotic value of /nl depends on the process by which the isocurvature modes 
become exhausted. We refer to this exhausted state as the adiabatic limit. Our arguments 
are phrased in terms of three-point correlations, but many of our conclusions are general and 
apply to arbitrary n-point functions including the two-point function. 

Classification of models. We restrict attention to models where a significant /nl is 
generated dynamically by inflation. This excludes examples such as the curvaton [28, 29] 
or modulated reheating [30, 31] which rely on an inflationary seed perturbation but amplify 
it by a noninflationary mechanism. If we disallow noncanonical kinetic terms, Maldacena's 
result guarantees that the models of interest must include two or more dynamically relevant 
fields [2]. 

Under these conditions, field fluctuations are generated at horizon crossing with almost 
Gaussian statistics [32-34]. However, unlike single-field inflation, these perturbations may 
cause spatially separated regions of the universe to experience different expansion histories. 
The set of phase space trajectories associated with an ensemble of such regions is initially a 
narrowly collimated bundle whose 'width' is set by quantum scatter. (We give details in §2.) 
The curvature perturbation is a precise measure of the relative expansion between spatial 
patches. Therefore its evolution is a consequence of focusing or defocusing of the bundle: 
as nearby patches of the universe evolve towards or away from each other in phase space, 
they experience varying expansion rates. In the adiabatic limit, the bundle degenerates to a 
caustic. Its width shrinks to zero, and all trajectories converge to a single line. 

When does convergence occur? The answer is model-dependent, but we can recognize 
broad classes of behaviour. 

• The potential may contain a focusing region, which enables trajectories to converge 
'naturally' If inflation ends in the vicinity of this region, an adiabatic limit is automatic. 
Examples include Nflation and related models [35-41]. However, one is always free to 
build models in which the adiabatic limit is achieved differently, perhaps by a waterfall 
transition. If the 'natural' limit is employed, there are two relevant questions. 

First, is the limit achieved before the end of inflation? If so, it is not necessary to 
specify details of the subsequent phases. Otherwise, we must choose among the various 
scenarios for reheating and later dynamics, and the predictions of the model may depend 
on our choice. 

Second, in the case where an adiabatic limit is reached during inflation, does this occur 
before the slow-roll approximation fails? This raises no issues of principle, but may 
influence how one chooses to study the model; for example, numerical methods may be 
required. We will return to this question in §5. 

• Alternatively, there may be no natural means by which trajectories converge. In such 
models there is no alternative: the inflationary model does not make unambiguous 
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predictions by itself, but only when embedded in a larger scenario which determines at 
least the mechanism by which inflation ends and the universe reheats. 

None of these observations are new, but their application to non-Gaussian statistics has yet 
to be studied in detail. 

Accepting the separate universe principle, to be discussed in §2 below, one obtains ex- 
plicit expressions for the re-point functions of the curvature perturbation [26] which automat- 
ically respect these conclusions. Therefore, concrete predictions can be obtained whenever it 
is possible to calculate these expressions accurately until the onset of conservation. However, 
in many cases this ideal procedure is impossible or impractical. Working in a special class of 
models where the potential is separable, Meyers & Sivanandam [42, 43] argued that the con- 
nected n-point correlation functions would be damped towards slow-roll suppressed values. 
Our analysis is closely related, but we argue that the value achieved in the adiabatic limit 
need not be especially small [41]. Indeed, in some cases, the adiabatic limit is associated 
with growth towards the asymptotic value, rather than decay. 

Objectives. In this paper we apply these ideas to the primordial bispectrum, and its 
amplitude /nl- There are three principal objectives. First, we illustrate that /nl can be 
sensitive to details of when and how the adiabatic limit is reached. Using examples, we 
demonstrate that — even in models where convergence occurs naturally — /nl may depend on 
the details of this process. Second, the calculations necessary to obtain a precise estimate of 
/nl can be technical, perhaps requiring recourse to numerical methods. For some models, 
simple techniques exist which allow a qualitative estimate of the evolution and asymptotic 
value of /nl- We outline these methods and apply them to simple examples. Third, we 
use numerical methods to perform a detailed study of the evolution of /nl in a selection of 
models. By themselves these calculations already reveal interesting patterns of behaviour, 
but also provide guidance regarding the asymptotics of models where an adiabatic limit is 
reached only through the intervention of post-inflationary dynamics. 

Outline. The plan of the paper is as follows. In §2 we discuss the separate universe ap- 
proach to perturbation theory in phase space, and the SN formalism. In §3 we discuss 
mechanisms for generating large evolving non-Gaussianity, and estimates for the maximum 
value and its dependence on initial conditions. §4 includes a brief account of analytic ex- 
pressions for estimating /nl, and discusses conditions necessary for this value to be large. §5 
reports a detailed numerical study of /nl in a selection of models with two or more fields. 
We conclude in §6. An Appendix contains details of some analytic calculations. 

2 Phase space description of slow-roll inflation 

Consider inflation driven by multiple canonical scalar fields 4>i (where i = 1, 2, . . . , iVf), self- 
interacting through a potential W(4>i). Defining Wj = dW/d(fti, the scalar equations of 
motion are 



14. In the "slow-roll" limit where e ^ 1 there is a dynamical attractor, allowing the decaying 



Inflation occurs when e = —H/H 2 
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mode to be discarded and restricting evolution to an iVf-dimensional submanifold II' on which 
(for example) the fa are determined in terms of the fa. The growing mode on II' satisfies 
3Hfa + W.i = 0. In the slow-roll limit it is possible to write e as a sum of independent 
contributions from each field, yielding e = Yli e i + 0(e? ), where the e, satisfy 

In simple models it may happen that the matrix rjij = MpWij /W also has small components 
\Vij\ « 1. 

Density fluctuations are generated by the inflationary background and can be measured 
by the curvature perturbation on uniform density spatial hyper surf aces, denoted £. For 
cosmological purposes its statistical properties are characterized by low-order correlation 
functions, of which the first two are 

(GdCka) = (2vr) 3 5(k 1 + k 2 )P(h), (2.3a) 
(C kl Ck 2 Ck 3 ) = (2vr) 3 5(k 1 + k 2 + k 3 )B(k u k 2 , k 3 ). (2.3b) 

The amplitude of the three-point function is usually measured in terms of a momentum- 
dependent parameter /nl 5 satisfying 

where i, j € {1, 2, 3}. The problem at hand is to calculate /nl- 

Bundles of trajectories. A comoving scale k is outside the horizon when k/aH < 1, 
where a is the scale factor and H = a/ a is the Hubble parameter. An overdot denotes a 
derivative with respect to time. During inflation H is approximately constant, whereas a 
is growing rapidly. Therefore k/aH rapidly becomes negligible a few e- folds after horizon 
crossing. Smoothing over a comoving scale somewhat larger than (aH)~~ 1 , widely separated 
spatial patches will locally evolve like an unperturbed universe up to small corrections. This 
is the separate universe picture [44-46]. 

An ensemble of smoothed regions picks out a collection of trajectories in phase space. 
If the ensemble is drawn from a spacetime region of finite comoving extent L then we can 
expect this collection to be narrowly collimated provided L is not too large. 1 We refer to 
this ensemble of clustered trajectories as a bundle. When making predictions for microwave 
background scales it may contain of order 10 6 trajectories or more [53]. In practice the 
analysis is simplified by working in a thermodynamic limit where the bundle formally contains 
an infinite number of trajectories. To avoid spurious infrared problems we should demand 
that these reheat almost surely in the same vacuum. 



1 Comoving quantities such as L are not physical, and can not appear in predictions for observable quantities 
[47]. The discussion in this paper is independent of L, but as a point of principle L should be removed from 
physical quantities by supplying a distribution function for the large-scale modes [48-50]. After doing so, all 
predictions depend only on physical scales. Similar conclusions have been obtained by a number of different 
methods [51, 52]. 
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The isocurvature modes label Fermi normal coordinates on II', adapted to the bundle. 
Each isocurvature field s has equation of motion s = and constitutes a conserved quantity 
[6, 54]. Together, these conserved quantities identify a trajectory. A particular location 
on each trajectory is specified by the integrated expansion iV ~ lna(t). This 'trajectory' 
approach can be traced to Hawking's formulation of perturbation theory [55], and was applied 
to inflation by several authors [44, 56-58]. An explicit description in terms of trajectories on 
IT' was given by Salopek [54] and Garci'a-Bellido & Wands [6]. Locally, N and the isocurvature 
fields generate a coordinate chart on II'. Fluctuations along the same trajectory generate 
the adiabatic mode, £ = 5N. Isocurvature modes represent fluctuations between trajectories. 
In an iVf-field model, the reduced phase space supports Nf — 1 isocurvature modes. 

Bundle sections. Consider a foliation of II' by submanifolds which are nowhere tangent 
to the bundle. An important example is foliation by surfaces of fixed energy density. We may 
use any such foliation to replace N as a label for length along the trajectories. Intersecting the 
bundle with an individual hypersurface generates a cross-section with coordinates inherited 
from the isocurvature modes. For example, working on uniform density hypersurfaces, the 
scalar fields take values §\. These "coordinates" are not independent but are subject to the 
constraint dW{4>1) = 0, leaving the expected 2Vf — 1 isocurvature labels. 

If the bundle has degenerated to a caustic then each hypersurface intersects the bundle at 
a unique point <$>\. Making a small change of trajectory 8(j)j earlier in the bundle's evolution, 
one will observe no change in (j>%. Therefore 5(j)^ = and we conclude d(j)^/d<j)* = 0. One can 
regard this as a requirement that physical predictions computed from statistical properties 
of the bundle become independent of its initial conditions. Speaking loosely, we describe this 
behaviour as an "attractor." The attracting trajectory is the caustic, and we will sometimes 
refer to it as the limiting trajectory. An ensemble of smoothed patches traversing the limiting 
trajectory differ only by their relative position within it, making £ conserved. 

This argument identifies regions where dcpf/d^j — > with regions where C is conserved. 2 
We describe it as the adiabatic limit. In what follows we will see that estimates of the decay 
rate of d^/drf)*- play an important role in analysing the adiabatic limit. As an example, 
consider the purely scalar dynamics associated with slow-roll inflation. A common type 
of limiting trajectory lies on a valley floor in the landscape generated by the potential. 
Descending into the valley, the mass-squared matrix associated with perturbations orthogonal 
to the direction of motion should be nondegenerate, with all eigenvalues large and positive. 
Taking the smallest eigenvalue of magnitude ~ mj_ one will generically find d(j)f/d(p* ~ 
e -m 2 ± N/3H 2 _ ^y e w jii discuss this example more carefully in §4. 

Bundle averages. These principles give a procedure to determine the statistics of £ in 
an adiabatic limit. One uses the bundle of trajectories to determine any required n-point 
correlation functions, and then imposes the requirement d<j)^/d<j)* — > 0. Unfortunately there 
is no unique way to do so. Different approaches to the limit correspond to different focusing 
mechanisms. As we have explained, focusing regions may occur naturally in some models; 

2 The condition that dtfif/dcfij — > is sufficient to show that that £ becomes conserved within the separate 
universe picture, but this argument does not demonstrate that it is necessary. 
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the examples studied by Meyers & Sivanandam [42, 43] are of this type. But whether or not 
a model naturally contains an adiabatic limit we may usually elect to impose a different one, 
perhaps by enlarging the field content to include a waterfall transition. We will study some 
possibilities below. 

It is first necessary to obtain the relevant correlation functions. When the slow-roll 
attractor is operative, the e-foldings along each trajectory can be expressed as a function of 
its initial conditions. Therefore N = N ((/>*). Expanding in the neighbourhood of a fiducial 
trajectory yields 

(^5N = Y, W + \Y. N ,iM S( f>j +•••> ( 2 - 5 ) 

i ij 

where Nj = dN/d(ft* and Njj = d 2 N / dcj)* dtp* . The 5(p* measure deviations from the fiducial 
trajectory, and will typically be of order the quantum scatter. After restriction to connected 
correlation functions there is no dependence on the arbitrary choice of fiducial trajectory. 

Eq. (2.5) enables the low-order correlation functions to be expressed in terms of the 
data iVj, N t ij, which can be computed in some models [37, 42, 59]. One finds (£C) = 
N iNj{8(j)i8<f)j)*i which determines the power spectrum (2.3a). Similarly, /nl can be written 3 
[26] ' 

jnl — c — ; 72 — • 



3 Transitory behaviour of /] 



NL 



It was explained in §§1-2 that our interest lies in adiabatic regions where all isocurvature 
modes are exhausted, preventing further evolution of £. This limit need not be achieved 
smoothly. For example, in hybrid scenarios the inflationary phase is suddenly destabilized by 
a waterfall transition, leading to abrupt convergence of trajectories. Although a convincing 
demonstration has not yet been given, in some circumstances the subsequent dynamics may 
preserve the value of /nl at the waterfall. This strategy has been invoked by various authors 
[61-64]. Some fine-tuning would be required to arrange |/nl| S> 1 at the transition point. 
Nevertheless, in these scenarios and others it may be misleading to focus exclusively on 
regions where phase space trajectories naturally converge. For this reason we pause to study 
the qualitative evolution of /nl ; whether or not we are close to an adiabatic region. We focus 
on scenarios where its value changes rapidly, before returning to focusing regions in §4. 

Under which circumstances should we expect the moments of the distribution function 
to change significantly? The distribution function describes how trajectories cluster around 
the core of the bundle. It is conserved under linear evolution, but is sheared or distorted 
on curved paths [65, 66]. These effects reshape the distribution function: even when it is 
initially Gaussian we expect probability to be relocated from the core to the outer layers of 
the bundle. This is associated with the generation of significant third- or higher n th -order 
moments. 



3 We are neglecting intrinsic non-Gaussianities of the dtp*. For slow-roll inflation with canonical kinetic 
terms these are negligible whenever /nl is large enough to be observable [37, 60]. 
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Curved paths can be generated by many choices of microphysics. We study only an 
especially simple mechanism. Where the potential's topography includes a ridge or valley we 
will encounter curved trajectories diverging from the ridge or converging into the valley floor. 
Examples of such trajectories have been studied by a number of authors [38, 63, 64, 67]. 

3.1 Ridges: Diverging trajectories 

We restrict attention to two-field models, which already capture the principal dynamical 
features, and label II' by coordinates and x- We assume a "ridge" or separatrix at x — 0. 
In the neighbourhood of an arbitrary point (<f>o,0) on the ridge the potential will generically 
have the form W ~ Wq + go(<p — 4>o)~ \ m \x 2 - The mass-squared m x is positive, and omitted 
terms are higher-order in <j) — 0o and x- These become relevant at some point after the 
trajectory has been ejected from the vicinity of x = 0. The trajectory x = is classically 
stable, although depopulated by quantum fluctuations [68-72]. 

Trajectories. Measuring length along each trajectory by the energy density, the evolution 
equations are 

1 dcj) g 



3M|d(iP) gl + {mlxf 
1 d X m 2 x x 



(3.1a) 



(3.1b) 



3M2d(#2) gl + imlxf 
According to (3.1b), a trajectory emanating from (H*,x*) and evolving to (H c ,x c ) satisfies 

^(Xl ~ xl) + 4 ln - = 3 ^p(^ 2 " H D- (3-2) 
2 m l x x* 

If \Xc\ IS M/ m x then the logarithm dominates and the trajectories disperse linearly in the 
sense Xc = X*D, where the growth factor D satisfies D = qP( h *~ h c) an d ft = 3(Mpm x /go) 2 . 
Nonlinear dispersion occurs in the region \x c \ ^ Is'ol/ 771 -^ where the quadratic term dominates. 
The transition between the two is the "turn," beyond which each trajectory is ejected from 
the ridge and rapidly evolves to \xd 3> bol/ 77 ^' ^ ^ ne ^ urn we nave 

[Xturnl ~ ^1, (3.3) 

mi, 



turn | 



which makes the kinetic energy in each field roughly equal, |</>turn| ~ |xt« 

This leads to the following physical picture. Trajectories which are still close to the 
ridge preserve their initial Gaussian profile. Trajectories populating the downhill-edge of the 
bundle quickly slide away, generating a heavy tail at large |x c |. In this region kinetic energy 
has greater relative importance, slowing the expansion rate and enhancing the frequency of 
excursions to large negative 5N . Therefore this mechanism will generate negative /nl from 
a Gaussian distribution. 

Whether a large negative amplitude is achieved in practice depends on the initial distri- 
bution of trajectories within the bundle, the nonlinear relation between the fields and £, and 
for how long the mechanism operates. Sufficiently far down the ridge the trajectories depend 
on the completion of W . Therefore, the approach to an adiabatic limit cannot be described 
by the techniques of this section. 
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8N analysis. We now translate to £ and repeat the analysis in the language of the 5N 
method [26]. Consider two trajectories originating well before the critical turning point, 
but initially separated by a distance (5(j)*,6x*)- It is useful to define 5 = m x \x/go\, where 
5* <C 1 indicates the initial point is very close to the ridge. In this region surfaces of constant 
energy density in II' practically coincide with surfaces of constant eft. Therefore, to bring this 
pair of trajectories to a common energy density H = H* requires a small excess expansion 
6N « (2ep- 1 / 2 ^/M P . The subsequent expansion history, measured to a surface H = H c , 
can be written N = N(H C ; H*,x*)- 

In what follows we work without loss of generality on the positive branch x > 0, and 
suppose (j) — 4>q and x remain sufficiently small that higher-order terms in the potential do 
not become relevant. Once this assumption fails, (p may acquire a nonnegligible dependence 
on x*> potentially invalidating our conclusions. Passing to the limit where <5</>* and 5x* 
become infinitesimal, we conclude that on arrival at H = H c the trajectories have experienced 
expansion histories which differ by 



dN 



3Hl 
9o 



+ ISrr&Ml d X * 



Hi 



H'i 



H 2 d(H 2 



bo + ( m lx) 



212 



X 



dx 
dx* 



n. 



(3.4) 



where the partial derivative is to be evaluated at constant H*. 

Invoking the chain rule, Eq. (3.4) determines all derivatives of N. We find 



N 



Hi 



H 2 d(H 2 ) 

212 



m bo + ( m lx) 



9l 



3(m2 x ) 2 / d X 



9o + {m 2 x) 2 \dx 



+ X 



H* 



d 2 X 
dx* 



H* 



(3.5) 



So far our considerations have been general. Prior to the turn, Eq. (3.2) makes d 2 x/dxt 
negligible whereas dx/dx* ~ x/x* is exponentially growing. In this region our assumptions 
make the integrands of Eqs. (3.4) and (3.5) positive, and therefore both iV x and N iXX are 
negative and decreasing. 

If m x is not too small, the integrals of (3.4) and (3.5) are dominated by their upper 
limits — where the exponential growth is maximized. Taking the initial evolution in x to be 
almost negligible, this requires 



m 2 x » - 



(H c /H*) 2 



(3.6) 



where the approximate equality applies if H c is at least a little smaller than H*. If x is to 
be sufficiently light to acquire a quantum fluctuation then m x <C H* , and if both conditions 
are to be compatible we must require f # Cl, A short calculation yields 



N 



?>m 2 x B. 2 c 



N 



■v 



XX 



9l 



A* 



X* 



(3.7a) 
(3.7b) 



This relation between N jX and N^ xx is a consequence of the exponential growth of x prior to 
the turn. 
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Initially, N >x and N xx are small in comparison with and N^p. In addition, N^ x ~ 
5*/Mp is constant and can safely be neglected. Therefore £ is dominated by the fluctuation 
in <fi, which is practically Gaussian. Using (2.6), we find 



5 /NL 



2e* + 



N 
N, 



■x 



3 m 2 x 



1 /iV, x 



1 + Gvf 



(3. 



While |iV x | <C l-ZV^I, the first term dominates and (3.8) gives |/nl| ~ e* < 1. As the trajec- 
tory moves away from the ridge the x*-derivatives become increasingly important whereas 
the (/^-derivatives are constant. When \N tX \ and l-ZV^I are comparable, /nl is dominated by 
the second term in (3.8). In virtue of (3.6) and the initial condition 5* <C 1, this is much 
larger than e* and causes a spike in /nl- Estimating the peak to occur when N t< j, ~ — N, 
we find 



x- 



/nl 



peak 



5* 



-0.3e; 



1/2 M P 



IX* 



(3.9) 



where rj x ~ MpW^ xx /W is the standard r/-parameter associated with %. In this expression 
and similar ones below, including Eq. (3.18), the numerical prefactor is uncertain by an 0(1) 
quantity which depends on the precise balance between and N ;X at the peak. 

On approach to the spike, Eq. (3.8) predicts that /nl is negative and growing like 
(Xc/x*) e - Subsequently, \ continues to increase and \N^ X \ eventually dominates \N^\. In 
this region £ is composed almost entirely of the x fluctuation. The non-Gaussianity becomes 
practically independent of 5* and decays like {Xc/x*)~ 2 ■ [These estimates of the growth rate 
and decay rate are valid before the turn, where \c is growing exponentially as described 
below Eq. (3.2).] 

Dropping numerical factors of order unity and using (3.3) to estimate /nl when the 
fiducial trajectory passes the turn, we find 



/nl| 



"of * %L 

c 



(3.10) 



If Eq. (3.10) is not invalidated by higher-order terms in the potential, it implies that the 
height of the spike is adjustable independently of /nl on entry or exit. Since the peak /nl 
occurs before most trajectories in the bundle reach the turning point (3.3), our analysis will 
apply provided these higher-order terms become relevant only after the turn. 

Scaling relations. Eqs. (3.8) and (3.9) give interesting scaling relations for the peak |/nl|> 
and for its growth and decay near the spike. Eq. (3.9) suggests that the maximum |/nl| 
attained during the spike has a practically universal power-law scaling for any potential 
which can be approximated by the coefficients go and m x up to the turn of the trajectories: 
for such potentials we should expect |/nl| oc with exponent v « 1. In §5 we will use 

numerical methods to study models which exhibit this scaling behaviour. 

3.2 Valleys: Converging trajectories 

The converse process occurs when a trajectory approaches a valley, where a bundle of trajecto- 
ries is nonlinearly focused rather than defocused. As above, we specialize to a two-dimensional 
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field space labelled by coordinates ((/), x) and suppose there exists a valley aligned with the 
X direction. In the neighbourhood of the valley we write W « Wo + W$ + W x , where Wo is 
a constant and 

= (3.11a) 

W x = goX + \m 2 x x 2 - (3.11b) 

If > m x then the slopes will be relatively steep in comparison with the valley floor. 
Omitted terms are higher order in and x> but become increasingly irrelevant as 0, x ~ * 0. 
Sufficiently far from = the motion is almost entirely in the direction. 

During descent into the valley, trajectories populating the uphill edge of the bundle 
experience a larger velocity in the orthogonal x direction compared to those lower down the 
slope. Therefore the uphill edge is compressed towards the centroid, generating a nonlinear 
distribution. The tail of the distribution is again on the downhill side, but in this case the 
tail has lower kinetic energy and enhances the frequency of excursions to a large positive SN. 
Therefore this mechanism generates a positive /nl- 



5N analysis. The evolution equations are 

1 d0 WL 



(3.12a) 



3M*d(H 2 ) (W£) 3 + (W£) 2 ' 

x (3.12b) 



In analogy with the ridge case, it is helpful to define a dimensionless measure of distance, 5, 
from the valley floor. We choose 5 = W^/W x , which measures the relative partition of kinetic 
energy between the fields. Our analysis applies when the trajectories begin from an initial 
position sufficiently high above the valley, where 5* 3> 1 and only the 0-field is in motion. In 
this regime, Eqs. (3.12a)-(3.12b) can be integrated to find 

ml 4>l (3 ' 13a) 

^ + 4=^ + 4)^)^ (3.13b) 

up to corrections of relative magnitude 1/5 2 . Eqs. (3.13a)-(3.13b) cease to be a good approx- 
imation no later than <5 ~ 1, when eft ~ x an d the kinetic energy in each field is approximately 
equal. In typical models this occurs at the turn. 

For J > 1, surfaces of constant energy density are practically aligned with surfaces of 
constant <f>. Adopting the methods of §3.1, we bring a pair of nearby trajectories separated 
by the displacement (50* , o~x*) to a common value of H, and write the number of e- foldings 
to a subsequent surface of constant energy density H c as N = N(H C ; H*,x*)- Passing to 
the limit of infinitesimal 50* and 5x*, and using Eqs. (3.13a)-(3.13b), we conclude that on 
arrival at H = H c , the trajectories have experienced expansion histories which differ by 

.„ Ml Mli,, f H ° H 2 d(H 2 ) f f^y\ 
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where 4> is to be understood as a function of H and we have introduced the mass ratio 
fi = m^/m^ < 1. Corrections to Eq. (3.14) are suppressed by 1/5 2 . 

Eq. (3.14) reproduces many features of the ridge analysis. The derivative N $ is constant, 
whereas |iV x | is initially zero but growing. Performing the integral, we find 



.,2 



where the "growth factor" $(x) satisfies 



Hx) = _ {nx + (x m _ i) + Wb_ (ar i _ i) + _ 1} . (316) 

We have assumed Wo dominates W*, but the generalization to other cases is straightforward. 
At x = 1 we have <1>(1) = 0. For x < 1 the dominant growing term depends on microphysical 
details of the model. Under our assumption Wq 3> W^* the dominant growth is initially from 
x . Inflation will not end naturally in a model of this type, so some other exit mechanism 
must be invoked. We will see examples of this kind in §5. On the other hand, if Wq/Wa,* ^ 1 
the logarithm will initially dominate. In either case, the asymptotic growth in the limit x <C 1 
is from x . Therefore, for a typical model <&(x) is a complicated function determined by a 
competition for dominance between the various terms. However, remarkably, in many cases 
the behaviour of /nl is almost independent of these complicated microscopic details. 
Differentiating (3.15), we find 

2 . /</> 2 \l N x 



N **=r' + T, A {Wit' (3 - i7) 

where A(x) = (x — 1) dln<I>(x)/dx. Eq. (3.16) implies that A is growing as <p decreases to 
zero. While \N tX \ is increasing towards \N^\ we find /nl increases, achieving a maximum 
value when |iV x | ~ \N^\. First, suppose the A-dependent term is subdominant at this time, 
which implies A < fj,5 2 . We find 

1/2 m^Mp 

/NLLeak ~ Vx* 6 * ~ °- 3e * ~ 21 \' ( 3 - 18 ) 

9 + rn x \X*\ 

which is independent of the growth rate (3.16) and the mass ratio fi. In this sense, the 
maximum value (3.18) is a "universal" phenomenon. When \N tX \ > \N<f)\ we find that /nl 
decays like at least until A ~ /x<5 2 , when it may stabilize as we will explain below. 

In analogy with the ridge, this sequence of growth and decay gives rise to a spike in /nl- 
Ultimately will decrease until 8 ~ 1, and the subsequent behaviour of /nl must be 

determined by different methods, such as those described in §4. Written in terms of the 
dimensionless measure 5, Eq. (3.18) coincides with (3.9) with the identification <5 va n cy = 
I/Bridge ■ In this language, the sign of /nl is inherited from the sign of rj x . 

Second, consider the "nonuniversal" case where the A-dependent term dominates (3.17). 
In this case, /nl increases until |iV x | ~ achieving a value larger than (3.18). Its precise 

value is set by the ratio A/5*, and may depend on details of the potential, including the mass 
ratio. When |iV x | > \N )t p\ the time dependence of /nl is set by A/<J>. Its precise scaling 
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depends on the dominant term in <5. In particular, if <3? ~ x" 1 then A/$ is approximately 
constant and /nl does not decay. In such cases, /nl exhibits a plateau and it may no longer 
make sense to speak of a spike at all. After the turn is completed, the nonlinear deformation 
of the bundle will partially relax, leading to decay of |./V XX |. The precise details, including the 
decay rate, are model-dependent. Eventually the fields reach equipartition of kinetic energy 
and this analysis breaks down. 

4 Asymptotic behaviour of /nl 

Whether a large |/nl| can be generated during an epoch of slow-roll inflation — perhaps from 
the "spike" mechanisms described above — is irrelevant unless it can be preserved in some 
adiabatic limit. The methods of §3 are insufficient to resolve this question. 

The potential may be such that a focusing region is naturally available. If inflation 
terminates in this region then the transitory evolution of /nl studied in §3 has no necessary 
connection with its final asymptotic value. In certain circumstances, where the focusing 
region can be analysed in detail, a relatively simple statement is possible. These are the 
scenarios studied by Meyers &: Sivanandam [42, 43]. One might have thought that the final 
/nl would depend only on the local shape of the potential in the focusing region, which 
will typically be a stable parabolic minimum. If so, the asymptotic value of /nl would be 
universal among all potentials sharing a similarly-shaped minimum. However, this is not the 
case. As we will explain, the asymptotic value of /nl generally depends on properties of the 
potential far from the focusing region. 

If multiple focusing regions are available, one must be selected by a combination of 
dynamics and initial conditions. To determine which possibility should be expected by late- 
time observers who map the anisotropy of the CMB requires an understanding of the infrared 
structure of the entire inflating volume [48]. This difficult "measure problem" remains un- 
solved. 

If a natural focusing region is not available, or is not selected, then one must be imposed 
and the entire analysis becomes significantly more complicated. In this case the transitory 
evolution studied in §3 may become relevant. We will have little to say about this possibility, 
although we investigate some numerical cases in §5. 

Natural focusing. In this section, we study models where inflation ends in a region of the 
potential where the trajectories are naturally focused. Broadly speaking, two possibilities 
exist. 

• The asymptotic value of /nl generated during focusing may be unobservably small, 
erasing any transiently large non-Gaussianity generated by spikes or other features. 
This possibility was emphasized by Meyers & Sivanandam [42, 43], who worked with 
a particular class of separable iVf-field models to be discussed below. However, other 
possibilities exist. 

• It may be possible to make the focusing process itself generate a large /nl by suitable 
choice of W. An example of such a model was given by Kim et al. [41]. (Indeed, in this 
model, /nl grows sharply during approach to the adiabatic limit.) 
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In principle, the behaviour of /nl in a focusing region could be determined from (2.5)-(2.6) 
by imposing the limit <9</>?/<9(/>* — > 0. Unfortunately, it is not known how to compute the U SN 
coefficients" iVj and iVy for an arbitrary model. Therefore a systematic discussion of this 
limit must apparently await future analytic developments. 

Explicit expressions for the 5N cofhcients are known only in very restricted circum- 
stances. Formulae for quadratic potentials were discussed by Lyth &: Rodriguez [26], Lyth 
& Alabidi [36] and Alabidi [38]. Later, Vernizzi & Wands [37] and Battefeld & Easther [73] 
gave expressions for an arbitrary sum-separable potential. Taking W = Vi((fii), one finds 





where dcp^/dcj)* satisfies 



d<f>< 




(4.1) 



(4.2) 



A similar expression for a product-separable potential W = JT^ Vi{4>i) was obtained by Choi 
et al. [74]. Comparable results for a general class of sum- and product-type potentials were 
given by Wang [59] and are summarized in the Appendix. It is also possible to take the 
Hubble rate to be separable rather than the potential [75, 76]. 

Focusing in a valley. A typical example of a focusing region is a valley of the potential 
landscape, perhaps terminating in a local minimum. For Nf fields, there are at least Nf — 
1 heavy directions with masses greater than the Hubble rate. Quantum fluctuations are 
suppressed in these directions, which prevent the bundle from diffusing up the sides of the 
valley. The steep slopes cause exponential convergence, and rapidly focus the bundle to a 
line. 

In the neighbourhood of the valley floor, we assume it is possible to choose coordinates 
on field space for which the potential approximately separates 



W 



where ip labels distance along the valley floor — which may be a light direction — and the Nf — 1 
fields s a are stabilized with masses m a = V" > H. By a suitable choice of coordinates we can 
arrange that (s a ) = 0. To describe a complicated valley it may be necessary to glue several 
such regions together. Focusing on the particular region described by (4.3), we denote the 
field values on entry to its domain of validity (p and s a . These will be functions of the initial 
fields <f>*. This initial point could generically occur far from the valley, where (4.3) need not 
be a good approximation. 

The heavy fields s a evolve according to 3Hs a = — TH a s a . After TV e- foldings from the 
point of entry, one finds 

Sa = s a (4>*)e-fo N v a (N')dN'_ (44) 

The total number of e-folds available within the valley is model-dependent. In a long valley 
the focusing may practically go to completion, making s a effectively zero. Alternatively, if 
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the valley rapidly terminates in a local minimum there may be insufficient time to focus the 
bundle completely. 

The fields (f>k can be written as linear combinations of <p and the s a , giving = 
1W + @k s a- The 7^ and /3% are constants, which depend only on the choice of separable 
coordinates used in (4.3). They are independent of the entry point (tp,s a ), which implies 

Therefore d^/d^j behaves like a linear combination of derivatives ds^/dcj)*. 

The number of e-foldings, N c (cj)*), which occur between the entry point ((p, s a ) and the 
surface c will usually depend on the initial point </)*. Assuming N c ((p*) does not exhibit 
a dramatic sensitivity to these initial conditions, Eq. (4.4) shows that ds^/d^ will decay 
exponentially as the trajectory settles into the valley. Potentials may exist which violate this 
condition, but we believe it will be satisfied for a majority of trajectories which flow over 
reasonably smooth potential landscapes. Where it is satisfied, this estimate of the decay rate 
applies once a trajectory has been captured by the focusing region, no matter what form the 
potential takes globally. 

In Eq. (4.5) the term will typically decay exponentially, because ~ s 2 a whereas 
decays less rapidly. Therefore Eq. (4.4) implies the derivatives d^/dcj)* decay at least 
as fast as the lightest isocurvature field. We conclude 4 

b l ~ e -f N Vs(N') dN> _ g-^JV ( 4g ) 



where rj s = mm{r] a } and N is the same quantity occuring in Eq. (4.4). The final equality 
applies if rj s is approximately constant during the focusing process. 

Separable potentials. In a globally sum-separable model, for which JVj satisfies (4.1), it 
may happen that (4.6) is sufficiently powerful to make the final "c-term" irrelevant. In these 
circumstances the correlation functions of £, including the spectrum and bispectrum, can be 
determined from the remaining term of (4.1), which depends only on boundary data at the 
initial time. For correlation functions among fields carrying comparable momenta of order k 
this is often taken to be the horizon-crossing time \rj\ ~ 1/fe, where rj is the conformal time. 
For this reason, the scheme has sometimes been called the horizon- crossing approximation 
[39-41]. Despite the name, we caution that this approximation does not consist in assuming 
that the perturbations are constant after horizon-crossing, but rather that their values in the 
adiabatic limit can be determined in terms of the shape of the potential there. A similar 
procedure can be applied in product-separable cases. 

It is less straightforward to estimate the minimum number of e-folds required to make 
the c-terms of (4.1) negligible. Although Eq. (4.6) gives information concerning the decay 
rate, the number of e-folds required to damp any contribution from the c-terms depends on 
their amplitude on entry to the valley. This is a function of each species' relative contribution 
to the energy density of the universe on the initial and final slices c and *, from which it 



The asymptotic notation x x y indicates that x and y share a common decay rate. 
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does not appear straightforward to draw general conclusions. However, since the isocurvature 
masses should be comfortably heavier than the Hubble scale, the parameter rj s will typically 
be much larger than unity. In these circumstances, rather less than O(10) e-foldings are 
usually required to accumulate a very substantial suppression of the c-terms. 

In the language of Meyers & Sivanandam [42, 43] , this damping of the c-terms is precisely 
the exponential suppression which they suggested would drive the bi- and trispectrum to slow- 
roll suppressed values. In the language of §2 it represents focusing of the bundle to a caustic. 
Our analyses are entirely consistent, but it is helpful to recall that the *-term in (4.1) need 
not be especially small. We briefly comment on this possibility at the end of this section. If 
that is the case, suppression of the c-terms can cause the correlation functions to experience 
a short phase of exponential growth as they approach their asymptotic values. Note that 
all these conclusions depend on the existence of a globally separable potential. We are not 
aware of a systematic study of the asymptotics of /nl in more general cases. 

One might harbour some reservations that the c-terms do not decay if the fields settle 
into a stable minimum, for which V^/Vl diverges. Near an arbitrary point, which can be 
chosen as the origin without loss of generality, Vk can be written Vk ~ A + B(f>k and Vk/V k ' 
approaches a constant. Near a minimum, one finds instead Vk ~ A + Bcjyf. Therefore 

V> HQ ~ 2B ( 7fc ^ + E a /«) V v * K) ' 

If A = the prefactor decays. Since the potential is sum-separable, we may always redefine 
all but one Vk to satisfy this condition. However, if the potential is not zero at the minimum 
then the remaining Vk must have nonzero A. Eq. (4.7) shows that we should choose the field 
4>k to have nonzero overlap with the direction of the valley floor, ie., 7^ 7^ 0. Under these 
circumstances the right-hand side of (4.7) still decays (although perhaps at a reduced rate), 
because by assumption <p decays strictly more slowly than any isocurvature mode. 

Large non-Gaussianity after natural focusing. Is it possible to obtain large /nl at the 
adiabatic limit? Working in a sum-separable potential, the foregoing discussion implies that 
the c-dependent terms in (4.1) may be discarded provided enough focusing can be achieved 
before the end of inflation. In general, the conditions required to achieve large /nl may 
still be complicated. However, a relatively simple picture emerges if we assume that iVj is 
large for one field (f> (or at most a few such fields) [41]. Therefore, V^/V^ at horizon crossing 
dominates the analogous terms for all other fields and /nl can be written 



5 V" 

/nl--| M ^ 4- 

6 Ki 



(4.8) 



In a single-field model the quantity V'ljV^ would be the inflationary 77-parameter. But in 
an assisted inflation the total potential may be much larger than Va, [77, 78]. Therefore rj^ 
can remain small, making (f> light at horizon crossing and causing it to acquire a quantum 
fluctuation by the usual mechanism, while V£ /V<j> can be appreciable. We study an example 
of this type in §5.1.2. In such models the sign of /nl is inherited from an "enhanced" r\- 
parameter, as in Eqs. (3.9) and (3.18), but unlike these cases the enhancement is measured 
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by the initial share of the energy density contributed by (f>, rather than the parameters <5 r id ge j 
^valley in Eq. (4.8) this enhancement factor is W/V^ S> 1. 

If several fields have comparable iVj, their perturbations contribute equally to £ at the 
adiabatic limit and dilute any non-Gaussianity by the same interference effect which leads to 
the central limit theorem. Therefore the largest values of /nl will be achieved where a single 
field has a dominant iVj. 

A similar discussion can be given for product-separable potentials. In this case the 
formulas depend solely on one field, labeled which must be the field still evolving at the 
adiabatic limit. Therefore /nl ~ 2e£ — r^ fc , and a large /nl would require a violation of 
slow-roll. We conclude that large |/nl| is not possible at the natural adiabatic limit in this 
class of models. 

5 Models 

The results of §4 show that, even for models where an adiabatic limit can be approached 
analytically, numerical calculations may be necessary to determine the degree of focusing 
which occurs near the end of inflation. In other cases there is simply no alternative. 

In this section, we report the results of numerical simulations and compare the outcome 
to the analytic theory developed in §§3-4. In appropriate circumstances we show that the 
simplified description of "spikes" obtained in §3 is an accurate match for full numerical 
simulations. We give examples where the focusing described in §4 goes to completion — 
making the "horizon crossing approximation" highly accurate — and others where it does not. 
A case of special interest occurs when the bundle would focus only slightly after the end of 
inflation. One might expect that the error in analytic predictions based on the strict adiabatic 
limit would be small, but it transpires that /nl can be rather sensitive to the details of the 
model. In models where there is no natural adiabatic region, reheating must occur before an 
adiabatic limit is reached. In these cases we perform a qualitative study of the dependence 
of /nl on the details of the reheating phase. 

Evolution after slow-roll using 5N. The phase space description of inflationary trajec- 
tories was discussed in §2. During multiple-field slow-roll inflation, each trajectory lies on a 
submanifold II' of the full phase space. 

In a model more general than multiple-field slow-roll inflation, extra coordinates will 
typically be required. First, if the slow-roll approximation fails then one must work on the full 
phase space IT rather than the attractive submanifold II'. Therefore new isocurvature modes 
are typically required to label the conjugate momenta 7Tj ~ fa. Second, matter species other 
than scalar fields may be included, perhaps to describe a phase of reheating. In such cases, 
the full phase space splits into a Cartesian product constructed from the phase space for each 
species, and suitable isocurvature modes labelling all these coordinates will be required. For 
thermalized radiation, a common choice is the temperature, T. 

Several numerical approaches exist to compute the statistics of the density fluctuation 
[65-67]. Here we take the simple approach of calculating the derivatives of N using a fi- 
nite difference scheme. This requires the slow-roll approximation at horizon crossing, where 



- 16 - 



initial conditions are set, but not subsequently. We have verified that our results are in- 
sensitive to changes in the step size of the finite difference scheme. Although slower than 
other approaches [65, 66], direct 6 N has the advantage of straightforward comparison with 
our analytic methods. Moreover, it requires only the evolution of an unperturbed universe, 
making a simple description of reheating — assuming thermal equilibrium and a single radi- 
ation fluid — easy to implement. These assumptions are at best quasi-realistic, but serve to 
indicate a plausible phenomenology. 

5.1 Two-field models 

5.1.1 Transitory models with interruption 

All two-field models documented in the literature which produce large, transient /nl exploit 
the spikes described in §3 [38, 63, 64]. Eqs. (3.8) and (3.9) show that to obtain large |/nl| 
from a ridge, one must tune the initial conditions so that x*/M~P ^ 1- Also, if this large 
/nl is to be preserved in the adiabatic limit, Eq. (3.10) implies that some mechanism must 
operate to end inflation before the majority of trajectories in the bundle encounter the turn. 
In two-field models with separable potentials, the parameter combinations required to ensure 
these conditions were given by Byrnes et al. [63]. The observables predicted in such models 
depend strongly on the choice of exit mechanism. In certain cases, such as two-field hybrid 
inflation, it is possible that /nl is not erased. In other cases this outcome appears unlikely. 

Two-field hybrid inflation. This model was studied by Alabidi & Lyth [38] and later by 
Byrnes et al. [63, 64]. The potential is 

W = \m\<? + '-mix 2 + \ (flj A 3 + sJxV) + (a 2 - v 2 f (5.1) 

where <j> and x are slowly-rolling fields, and a is a waterfall field which becomes destabilized 
when g^cft 2 + g x 2 X 2 = Au 2 . We take the masses and m x to be positive, and assume 
9<t>/dx = m2 s/ m x- This ensures that the waterfall occurs at fixed energy density, making it 
unnecessary to account for the effect of inflation ending on different hypersurfaces [61, 62, 79]. 
If the masses are not equal, there is a steep slope in the direction of the more massive field. 
The trajectories evolve along this steep direction and then turn towards the global minimum. 
We expect some non-Gaussianity to be generated during this process. 

In Fig. 1 we plot the evolution of /nl for m ( j ) /m x = 5, 7)$ = 4Mpm 2 /(A^ 4 ) = 0.08 and 
initial conditions x* = O.OOIMp, 0* = 0.5Mp. We adjust the remaining parameters so that 
the waterfall occurs when |/nl| > 1 and takes much less than a Hubble time to complete. The 
blue dotted line represents the /nl generated by the slow-roll fields, ignoring the waterfall. 
If inflation fails to terminate on the spike this leads to a negligible asymptotic bispectrum. 
The potential is of 'valley' type with g = 0, for which the analysis of §3 explains the positive 
spike. Beginning from /nl ~ e* there is rapid growth to a positive peak, followed by a softer 
decay. The peak value is well approximated by Eq. (3.18) which yields /nl ~ 8. By varying 
the initial conditions, we have confirmed that the peak value scales approximately as 1/x*, 
as predicted by (3.18). 

The solid red line represents a numerical evolution, terminated by a waterfall transition 
on the growing arm of the spike. In the early stages, the numerical results follow the analytic 
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Figure 1. Evolution of /nl for the model (5.1) calculated numerically (solid red line) with the hybrid 
transition included for the parameter values in the text, the (red) dot-dashed line is added to illustrate 
the nearly constant final level of /nl- The (blue) dashed line, represents the analytical evolution, with 
no hybrid transition included. 



prediction. We continue the calculation past the end of inflation by allowing the waterfall field 
a to become operative. Care must be taken in modelling this transition. We give the waterfall 
field a small value consistent with the typical RMS value expected from quantum mechanical 
excitations of a massive field in de Sitter, crms ~ H 3 /M, where M is representative of 
the waterfall mass before the transition. The results are extremely insensitive to its precise 
value. 5 

The hybrid field is heavy at horizon crossing and is therefore unperturbed. Hence, 
we need not differentiate iV with respect to a. Using these assumptions, we find that /nl 
appears to be conserved through the hybrid transition (see Fig. 1). The numerical evolution 
is only continued for a fraction of an e-fold after the transition, during which time /nl does 
evolve due to oscillation of the primary fields. However, the resulting oscillations in /nl are 
decaying and are centred around a fixed value. This behaviour is apparently generic for a 
range of parameter values, provided the transition happens sufficiently rapidly — in less than 
an e-fold. One must already impose this "rapid transition" condition to avoid issues with 
primordial black holes [80]. 6 

5 In reality, this RMS value is made up of many inhomogeneous short scale modes. Their collective evolution 
approximates that of a homogeneous mode, at least in the initial stages before the minimum is reached [80]. 
We expect this approximation captures at least some of the physics which occurs at the hybrid transition. 

6 Note that when the hybrid transition does not occur on a uniform density hypersurface a significant 
extra contribution to /nl can be generated. We have checked a small number of these cases numerically and 
find agreement with the formulae given in Ref. [64], although we have only considered examples of positive 
curvature rather than ridges. 
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5.1.2 Large non-Gaussianity at the natural adiabatic limit 

We illustrate this case using a model closely related to the iV-axion model of Kim et al. [41], 
in which the potential is taken to be V = Yli^fi^ ~ cos 2vr/~ 1 (/)j). The sum is taken over 
a large number of uncoupled axions, and fi is the decay constant for the i th axion. We will 
study further examples of this type in §5.2. In Ref. [41], many axions were invoked to generate 
a phase of assisted inflation. Because the potential is sum-separable, the perturbations can 
be calculated at the adiabatic limit using (4.1) after dropping the c-term provided inflation 
ends when the final field gracefully exits from slow-roll. Whether this occurs depends on the 
number of fields and the choice of fi. Taking fi = f< Mp for all i and supposing that the 
initial conditions are chosen so that only a small number of fields populate the hilltop region 
near fa = 0, the asymptotic /nl can be calculated using (4.8). It will typically be moderate 
or large. 

In this section we study a related two-field model. Dynamically, the large number of 
axions which begin away from the hilltop region serve only to source the Hubble rate. The 
single field closest to the hilltop sources the non-Gaussianity. (This model has some similarity 
to the scenario of Boubekeur & Lyth [81].) Therefore, most of the axions can be replaced by 
a single effective field with a quadratic potential, retaining the full cosine only for the axion 
closest to the hilltop, 



where A and / are constants. 

Near the hilltop, the axion potential approximately satisfies V(x) = 2A 4 (1 — tt 2 x 2 // 2 )- 
This yields a tachyonic mass 2ttA 2 / f. Adjusting the <f) potential if necessary to ensure that 
X remains light at horizon crossing, the mass induces a large /nl via (4-8). 

In Figs. 2 and 3 we show a numerical evolution for two choices of parameters. In Fig. 2 we 
take / = Mp and A 4 = 25m 2 / 2 /(47r 2 ), which makes the mass of the axion five times greater 
than the mass of <f>. The initial conditions are 0* = 16Mp and x* = (//2 — 0.001)Mp. As a 
consequence of its large mass, the axion rolls off the ridge quite early. Therefore the system 
evolves to the limiting trajectory long before the end of inflation. According to §3, departure 
from the ridge should produce a large negative spike. Later, convergence into the minimum 
should produce a positive spike, perhaps followed by a plateau. Finally, as the isocurvature 
modes are exhausted, the system should evolve to the adiabatic limit (4.8). These features are 
clearly visible in Fig. 2. Matching the potential (5.2) to the analysis of §3 and using Eq. (3.9), 
we expect the negative peak of /nl to occur at /nl ~ —0.3Mpel^ 2 (f /2 — x*)" 1 ~ —26. This 
gives good agreement with the observed value. By varying the intial conditions we have 
verified that scaling with l/(//2 — x*) _1 is reproduced to a good approximation. In this 
case the difference between our slow-roll analysis and the full numerical calculation is at the 
level of a few percent, consistent with the accuracy of the slow-roll approximation. 

In the second example we take / = Mp and A 4 = m 2 / 2 /(47r 2 ), giving both fields the 
same mass. In this case the axion starts to evolve only near the end of inflation, where <fi is 
approaching the minimum. Indeed, much of its evolution takes place while <j> is oscillating. 
In these circumstances the adiabatic limit cannot be calculated analytically using (4.1). 
Nevertheless, the important features can still be understood. While 4> is oscillating in the 




(5.2) 



- 19 - 



20 
10 


-20 
-30 

~ 4 °0 10 20 30 40 50 60 
N 

Figure 2. Evolution of /nl for the model of Eq. (5.2) (first set of parameter choices). The solid red 
line is a numerical calculation. The blue dashed line is an analytic prediction. The horizontal green 
dashed line represents the analytically calculated adiabatic limiting value. 

minimum, its potential energy contributes to the energy density in a way not accounted for by 
the slow-roll approximation. If we suppose the cf) oscillations do not lead to rapid reheating or 
preheating, we may expect £ to approach a constant as the trajectories settle in the minimim 
and Hubble friction drains their energy. The results are given in Fig. 3. In this simple 
example, /nl oscillates around an asymptotic value which is lower than would be expected 
if the adiabatic limit were reached during inflation. In more sophisticated examples, where 
complex dynamical behaviour can occur during the oscillating phase, it would be necessary 
to follow their decay in precise detail [82-88]. 

This example is representative of a class of model where natural focusing occurs — in this 
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Figure 3. Evolution of /nl for Eq. (5.2) (second set of parameter choices). The solid red line is 
a numerical calculation. The blue dashed line is an analytic prediction. Only the final few e-folds 
are shown. The horizontal green dashed line represents the analytically calculated adiabatic limiting 
value. The solid vertical line indicates when inflation ends, computed using the exact equations of 
motion. As the axion rolls, inflation momentarily restarts and the slow-roll expressions cease to be a 
good approximation. 





- 20 - 



case, caused simply by Hubble damping — but does so only after the slow-roll assumption is 
violated. There are other models in this class which lead to a large non-Gaussianity at the 
natural adiabatic limit, such as models which possess an inflection point in their potential 
with a very slight gradient. We intend to return to these cases in future work [89]. 

5.1.3 Models with no reconvergence in field space 

The third possibility discussed in §1 occurs when the trajectories disperse in field-space but 
the potential provides no region which would enable them to refocus. An example is provided 
by the model 

V = V^ 2 e-^\ (5.3) 

which was introduced by Byrnes et al. [63]. The dynamics were followed only until the end 
of slow-roll, at which time the non-Gaussianity was indeed large. However, at this point, the 
isocurvature modes were not exhausted and the curvature perturbation was still evolving. 
To study this model, we make the the same parameter choices as Byrnes et al., setting 
A = 0.05/M2, fa = 16M P and X i = 0.001M P [63]. 

The initial stage is descent from a ridge, and therefore we expect /nl to approach 
negative values. This is confirmed in Figs. 4 and 5. As the bundle rolls along the ridge 
(defined by x = 0), we find that Eq. (3.9) gives /nl ~ —26, in good agreement with the first 
negative peak in /nl- We have confirmed that the expected scaling with initial conditions 
is approximately respected. With this choice of parameters, a large /nl is still present 
as slow-roll breaks down. But because no limiting trajectory is available, /nl continues to 
evolve — and subsequently oscillates wildly, as the fields oscillate about the line (ft = 0. This 
does not represent a stable attractor: <ft = is a degenerate vacuum, and the field evolves only 
along the x direction during the oscillations. To reach an adiabatic limit we must apply a 
prescription for reheating. Here, we adopt a very simple perturbative model in which energy 
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Figure 4. Evolution of /nl for the model (5.3). The solid red line is a numerical calculation for the 
parameter values quoted in the text and T = 0. 
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Figure 5. Upper panel: evolution of /nl calculated numerically (solid red line) for the parameter 
values quoted and T — (Vo/10) 1 / 2 /Mp. The blue dashed line represents the corresponding plot 
with T = {Vo/IOOY^/Mp, and the green dot-dashed line represents the analytical evolution. The 
solid vertical line indicates the end of inflation, measured using the analytic evolution. The analytic 
evolution terminates when the x field reaches zero, as the slow-roll expressions can't evolve the fields 
past this point. Lower panel: magnified in the vicinity of the end of inflation. 



is transferred from the field into a radiation component. The dynamical equations are 

dW 

4>i + 3H(j>i = Tifa - — (5.4a) 
o<Pi 

p = -4Hp + J2^l (5.4b) 

i 

where p is energy density of radiation, and the Tj represent the decay rate from species i. 
We illustrate the effect of reheating in Fig. 5. The final value of /nl is sensitive to the choice 
of Tj, and hence the time-scale of reheating. We take Tj = T for all i, making reheating 
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begin approximately when H = T and take place on a uniform density hyper surf ace. A more 
complicated prescription leads to strong secondary effects which radically alter the value or 
sign of /nl- After reheating, if the radiation is the only contribution to the energy density, 
then the statistics of £ at this time will be the ones relevant for observation. 

Our aim has not been to present a realistic model. Rather, we wish to demonstrate 
that, if no attractor exists within the inflationary regime, we must follow the dynamics until 
all observable quantities stop evolving at the adiabatic limit. We can expect the asymptotic 
value of each observable to be sensitive to this evolution, including the time scale and details 
of reheating. 

5.2 iV f -field models 

Similar results naturally apply in models with a larger number of fields. In this section we 
study the model of Kim et al. [41] involving many axion fields self-interacting though the 
potentials 



Where the parameters Aj and /, take common values A and / for each species, and / < Mp, 
this generates naturally large /nl at the adiabatic limit. Although larger /nl can naively 
be obtained by decreasing the fi, it is necessary to simultaneously increase the number of 
fields in order to obtain sufficient inflation. There is another difficulty. As the fi decrease, 
the approach of /nl to its asymptotic limit occurs later in the evolution. In Fig. 6 we show 
one realization of this behaviour for Nf = 1800 and fi = Mp, with initial conditions for the 
fields randomly distributed in the range < fa < irMj>. The slow- roll phase ends at latest 
when e = 1, marked by the vertical black line. The evolution to the right of this line is not 
trustworthy and should be replaced by a numerical calculation. Unfortunately, owing to the 
large number of fields we have not been able to perform a non-slow roll analysis due to the 
prohibitive running time of the computation. 

There is a specific case where this model can be related exactly to the two-field axion- 
plus-quadratic model of the the previous section: when 2Vf — 1 fields are initially close to 
the minimum of the axion potential (fa <C //2 for i = 1, . . . , N{ — 1) with identical initial 
conditions, and one field is initially close to its maximum, <pN t ~ //2. In this case the Nf — 1 
fields act like a large number of fields with a quadratic potential. When they all evolve from 
an identical initial condition, the dynamics of the many fields is completely identical to the 
dynamics of a single field <I> 2 = X^i" 1 ^?' wr th a quadratic potential of the same mass as 
the individual fa fields. For / = Mp, this reduces identically to the second of the two-field 
axion-plus-quadratic cases studied in §5.1.2. (See Fig. 3.) 

Finally, it is interesting to note that the Af-field axion model — for which a large /nl 
follows from relatively generic initial conditions — is closely related to a two-field model in 
which generation of large /nl apparently requires significant fine-tuning. The tuning appears 
less dramatic in the original Af -field model. It is interesting to conjecture that the fine-tuning 
of initial conditions required to give large /nl in two-field models may be reduced in models 
with many fields. 




(5.5) 
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Figure 6. Evolution of /nl for the iVf-axion model, calculated analytically (solid blue line) under 
the slow-roll approximation. The horizontal dot-dashed red line is the asymptotic value computed 
using the horizon-crossing approximation. (This is unreliable in the present case.) The vertical black 
line corresponds to e = 1. Since /nl has not reached the adiabatic limit at this point, this model is 
an example in which the adiabatic limit is reached after slow-roll ends. Therefore, numerical analysis 
is required to obtain a reliable value for /nl- 

6 Conclusions 

We have studied the evolution of non-Gaussianity in multiple-field models of inflation. Unless 
all isocurvature modes become exhausted before the end of inflation, we find that there need 
not be a unique prediction for /nl- Instead, the final value can depend on independent details, 
such as the microphysics of a reheating or preheating phase. Where the trajectories naturally 
focus — for example, if inflation ends with all fields settling into a minimum of the potential — 
numerical calculations are typically required to determine the precise asymptotic value for 
/nl- If there is no natural focusing region then numerical calculations and a prescription for 
reheating will be required. This confirms the natural expectation that analytic predictions 
[36-38, 73] are reliable only if the flow of power from isocurvature to curvature modes is 
quenched before the end of the slow-roll phase. 

If an adiabatic limit is reached without passage through a natural focusing region, 
perhaps by invoking a waterfall transition, then this may occur when the value of /nl is 
transiently large. However, we caution that although our numerical calculations indicate 
that /nl can sometimes be preserved through a hybrid transition, there does not yet appear 
to be a precise characterization of the conditions required for this to occur. Also, whether 
the end-point of the waterfall is an adiabatic limit may be model dependent. Temporarily 
ignoring these subtleties, we have shown that descent from a ridge or convergence into a 
valley can result in a significant, dynamical but transient enhancement of /nl- The two 
cases are distinguished by a different sign of the resulting /nl> which is inherited from the 
local r\ parameter. 

We have verified that it is possible to construct sum-separable models which exhibit 



large /nl even when the adiabatic limit is reached during slow-roll inflation. Therefore there 
is no correlation between large /nl at the end of inflation and the presence of an inexhausted 
isocurvature perturbation, as has occasionally been suggested. On the other hand, we have 
demonstrated that this is impossible for product separable cases, where /nl is always of order 
the slow-roll parameters at horizon crossing. 

Among the models we have studied is a new two-field model related to the Af-field axion 
model [41]. This exhibits large /nl at the adiabatic limit when this limit is reached before 
the breakdown of slow-roll. In addition, inflation can end gracefully rather than through a 
sudden transition. As far as we are aware, this is the first example of such behaviour in the 
two-field context. This model explicitly illustrates that predictions of the A r f -axion model 
may be modified if a full numerical calculation for a sufficiently large number of fields could 
be performed. 
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A Appendix 



A.l Calculating JVj 

Under the assumption of slow-roll and monotonicity {fa < 0), the number of e- folds can be 
written with the field fa as a time variable. Taking the functional derivative of this integral 
generates three components — two "boundary terms" evaluated on the initial ('*') and final 
('c') slices, and a "path term", 

d 



N,iMl 



W 



Wj 



w 



W, 



w 



dfa. 



(A.l) 



The summation convention is not used. Physical quantities are independent of k, which 
may be chosen arbitrarily. Employing the notation S = ^ an d P = Yli=i Vii where 

V\ = V\{4>i) are dimensionless, we restrict attention to potentials of the form W = Mp F(S) 
and W = Mp G(P), where F and G are arbitrary functions. We refer to these as sum- and 
product-separable potentials respectively. 

For the product-separable potential W = Mp G(P) we evaluate the path term and the 
final boundary term in Eq. (A.l). Using the slow-roll parameters y/2ei = MpG'PV-/GVi 
and also defining m = e^/e we find 



_L Yl 

Ml VI 



Gm 



G'P 







G 



(A.2) 



This applies for arbitrary G. Note that any dependence on k has disappeared. If the ratio 
G/G'P depends on P, then the integral requires knowledge of the variation of Vj with fa. 
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There are at least two cases in which the P-dependence is lost: if the integrand is either zero 
or constant. If the integrand is zero then we have G/G'P = A for some constant A. This 
gives the general solution G = AP B , where the constants A and B can be absorbed into a 
redefinition of the potential. This yields 

1 vS 



N, 



W = Mp P. 



(A.3) 



If the integrand is constant this implies G = (B\nP + D) l l A , where A,B,D are constants 
of which B and D can be eliminated by a further redefinition. We conclude that a general 
potential yielding a P-independent integrand can be expressed as W = Mp (lnP) 1 ^ and 
gives 



_A_ Vi_ 
Ml V> 



(In V* - In V? + (In P) c <) , W = M 4 (In P) 



,1/A 



An identical procedure applies to the sum-separable potential, leading to 
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N, 



U; 



A 



W 



M*e s 



-(V*-Vf + S c ul), W = Ml,S l l A . 



(A.4) 

(A.5) 
(A.6) 



A. 2 Correspondence between different separable potentials 

There is a correspondence between Eqs. (A.5) and (A.3), and between Eqs. (A.6) and (A.4) 
which was first pointed out by Wang [59] for two field potentials. Redefining In V{ — > Vi turns 
a W = Mp P potential into a W = Mp e s potential, and redefining e Vi — > Vi transforms the 
potential W = M| S l/A into the form W = Mp (InP) 1 ^. Therefore, it is unnecessary to 
proceed with all four classes of potential. As an independent pair of potentials, we choose 
W 



Mp P and W 
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Mp-P: 
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For W 



M* S l ' A : 




uf) and A^., 
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MpS*^/2e 

These formulae have previously appeared in Refs. [59, 63, 73]. 
A.3 Nij and / NL 

To calculate higher-order statistics we require the second derivatives iVy 
nient to introduce new dimensionless slow-roll parameters 



(V*-Vf + S c ut). (A. 
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These are elements of a vector a and matrix b respectively. Note that b is diagonal for the 
product-separable potential W = Mp P and so we rewrite Ni = u-/Mp a* with U{ = af/\a\ 2 . 
A simplification occurs if the limiting trajectory is a straight line in field space, which implies 

/NL| straight = -qT^U ^Z^if^ ( A - U ) 

i=l 1 

where the vector q c has elements qf = 2Vj. If the limiting trajectory lies along the <j) axis, then 
/nl = — I n this case, the sign of /nl is given by the mass of the field at horizon crossing 
but the magnitude is slow-roll suppressed. If we do not make the simplifying assumption 
that the limiting trajectory is a straight line, then 



i N t / c3 \ 



a • b • a 

+ 2 , , n 



(A.12) 
(A.13) 



For the potential W = Mp S 1 ^ we have N { = (V* - Vf + S c ul) /M P a*S* . In the 
simpler situation where u\ is a constant, we find 

«P ~«l»» lg h, = ASa - qi (Aa- + J) - - „J) (A.14) 

/NL I straight = f - f ^ £(* (*< + f ) + " ^ 

Dropping the assumption that u\ is constant leads to additional terms of the form 

N ,ij = ^straight + M$-^ a * a *lac\2 ( A - 16 ) 

" /NL| straight + 3^2 | c|4| a c|2 -rf^" ^-17) 

i,j=l J 

Mij = dibijdj — aiUj(b • a)i — ajUi(b • a)j + UiUj(a • b • a). (A. 18) 
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